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Abstract 

We provide the first analytical results for the connectivity of dynamic random 
geometric graphs - a model of mobile wireless networks in which vertices move 
in random (and periodically updated) directions, and an edge exists between two 
vertices if their Euclidean distance is below a given threshold. We provide precise 
asymptotic results for the expected length of the connectivity and disconnectivity 
periods of the network. We believe the formal tools developed in this work could 
be of use in future studies in more concrete settings, in the same manner as the 
development of connectivity threshold for static random geometric graphs has 
affected a lot of research done on ad hoc networks. In the process of proving 
results for the dynamic case we also obtain asymptotically precise bounds for the 
probability of the existence of a component of fixed size £, £ > 2, for the static 
case. 



1 Introduction 

Random Geometric graphs (RGG) have been a very influential and well-studied model 
of large networks, such as sensor networks, where the network agents are represented 
by the vertices of the RGG, and the direct connectivity between agents is represented 
by the edges (see for example the recent books 0[12]). Informally, given a radius r, 
a random geometric graph results from placing a set of n vertices (agents) uniformly 
and independently at random on the unit torus and connecting two vertices if and 
only if their distance is at most r. 

In the late 90s, Penrose p3], Gupta-Kumar [7] and Apple and Russo [1], gave ac- 
curate estimations for the value of r at which with high probability, a RGG becomes 
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connected (see Section 2). We denote this value of r by r c . Thereafter, hundreds 
of researchers have used those basic results on connectivity to design algorithms for 
more efficient coverage, communication and energy savings in ad hoc networks, and 
in particular for sensor networks (see the previously mentioned books). On the other 
hand, much work has been done on the graph theoretic properties of static RGG, 
comprehensively summarized in the monograph of M. D. Penrose |15| . In Section [21 
we prove a result on static random geometric graphs, which was not known before 
(Theorem 1): At the threshold of connectivity r c and for any fixed I > 1, the proba- 
bility of having some component of size at least I other than the giant component is 
asymptotically 0(1/ log^ -1 n). Moreover, the most common of such components are 
cliques with exact size I. This result plays an important role in the derivation of the 
main result for the dynamic setting, which is explained below. 

Recently, there has been an increasing interest for MANETs (mobile ad hoc net- 
works). Several models of mobility have been proposed in the literature - for an 
excellent survey of those models we refer to [ID]. In all these models, the connec- 
tions in the network are created and destroyed as the vertices move closer together 
or further apart. In all previous work, the authors performed empirical studies 
on network topology and routing performance. The paper [5] also deals with the 
problem of maintaining connectivity of mobile vertices communicating by radio, but 
from an orthogonal perspective to the one in the present paper - it describes a kinetic 
data structure to maintain the connected components of the union of unit-radius disks 
moving in the plane. 

The particular mobility model we are using here (in the literature it is often called 
the Random Walk model) was introduced by Guerin [6], and it can be seen as the 
foundation for most of the mobility models developed afterwards [10J. In the Random 
Walk model, each vertex selects uniformly at random a direction (angle) in which 
to travel. The vertices select their velocities from a given distribution of velocities, 
and then each vertex moves in its selected direction at its selected velocity. After 
some randomly chosen period of time, each vertex halts, selects a new direction and 
velocity, and the process repeats. An experimental study of the connectivity of 
the resulting ad hoc network for different values of n and r for this particular model 
is presented in [17]. As is stated in the same paper, in many applications which are 
not life-critical connectivity is an important parameter: "Temporary network discon- 
nections can be tolerated, especially if this goes along with a significant decrease of 
energy consumption. " In the present paper, we perform the first analytical study 
of connectivity in the Random Walk model. The particular variant of the model, 
that we study, is the following: Given an initial RGG with n vertices and a radius 
set to be at the known connectivity threshold r c (see Section [2]) , each vertex chooses 
independently and uniformly at random an angle a E [0, 2ir), and moves a distance s 
in that direction for a period of m steps. Then a new angle is selected independently 
for each vertex, and the process repeats. We denote this graph model the Dynamic 
Random Geometric Graph. 

Our main result (Theorem 13. ip provides precise asymptotic results for the ex- 
pected number of steps that the graph remains connected once it becomes connected, 
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and the expected number of steps the graph remains disconnected once it becomes 
disconnected. These expressions are given as a function of n and s. Surprisingly, the 
final expression on the length of connectivity does not depend on the size of the inter- 
vals between changes of angles (as long as the angles do change, after some possibly 
large number of steps). It is worth to note here that the evolution of connectivity in 
the dynamical setting described in Section [3] is not Markovian, in the sense that stay- 
ing connected for a large number of steps does have an impact on the probability of 
being connected at the next step. However, one key (and rather counterintuitive) fact 
is that, despite this absence of the Markovian property, our results, by Lemma l3.6j, 
depend only on the probability of connectivity change in two consecutive steps. 

Throughout the paper, we assume the usual Euclidean distance on the 2-dimensional 
torus [0, l) 2 , but similar results can be obtained for any £ p -norm, 1 < p < oo. More- 
over, our results also can be extended to any cube [0, l) d for any d = 6(1). To the 
best of our knowledge, the present work is the first work in which the dynamic con- 
nectivity of RGG is studied theoretically. In [I] the loosely related problem of the 
connectivity of the ad hoc graph produced by w vertices moving randomly along the 
edges of a n x n grid is studied. In [11] the authors use a similar model to the one used 
in the present paper to prove that if the vertices are initially distributed uniformly at 
random, the distribution remains uniform at any time. Unless otherwise stated, all 
our stated results are asymptotic as n — > oo. As usual, the abbreviation a.a.s. stands 
for asymptotically almost surely, i.e. with probability 1 — o(l). All logarithms in this 
paper are natural logarithms. 

2 Static Properties 

In this section, we give the basic known results about connectivity of static RGG, 
and provide a new bound on the probability of the existence of components of size 
i, i > 2, at the connectivity threshold of a RGG. The threshold for connectivity of 
RGG has a long and exciting history, but due to lack of space, we only refer to [15]. 
The formal definition of random geometric graph is the following: Given a set of n 
vertices and a positive real r = r(n), each vertex is placed at some random position 
in the unit torus [0, l) 2 selected independently and uniformly at random (u.a.r.). We 
denote by X.- L = (xj,t/i) the random position of vertex i for i E {1, . . . , n}, and let 
X = X(n) = Ur=i^«- We note that with probability 1 no two vertices choose the 
same position and thus we restrict the attention to the case that \X\ = n. We define 
G(X; r) as the random graph having X as the vertex set, and with an edge connecting 
each pair of vertices Xi and Xj in X at distance d(Xi,Xj) < r, where d(-,-) denotes 
the Euclidean distance in the torus. 

Let K\ be the random variable counting the number of isolated vertices in G(X; r). 
Then, by multiplying the probability that one vertex is isolated by the number of 
vertices we obtain E (Ki) = nil — 7rr 2 ) n_1 = ne~ 7TT n — 0(r 4 n). 

2 

Defining fi = ne~ 7TV n , we get that the asymptotic behaviour of fi characterizes 
the connectivity of G(X;r): 
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Theorem 2.1. 



• IffA—^O, then a.a.s. G(X;r) is connected. 

• If fj, = ©(1), then a.a.s. G(X;r) consists of one giant component of size > n/2 
and a Poisson number (with parameter fi) of isolated vertices. 

• If fj, —* oo, then a.a.s. G(X;r) is disconnected. 

From the definition of /i we deduce that fi = ©(1) iff r = J ^^ ^. Therefore 
as a weaker consequence we conclude that the property of connectivity of G(X;r) 

exhibits a sharp threshold at r = \ f^j^ L - Theorem 12.11 also implies that, if /i = ©(1), 
the components of size 1 (i.e. isolated vertices) are predominant and have the main 
contribution to the connectivity of G(X;r). In fact if C (respectively T>) denotes the 
event that G(X;r) is connected (respectively disconnected), we have the following 

Corollary 2.1. Assume that /i = 0(1). Then 

P(C) ~ P(Ki = 0) ~ e-», P(D) ~ P^i > 0) ~ 1 - e't 1 . 

Therefore, if /i = 0(1), the probability that G(X;r) has some component of size 
greater than 1 other than the giant component is o(l). The goal of the remainder 
of this section is to give more precise bounds on this probability. We need a few 
definitions. 

Given a component T of G(X;r), T is embeddable if it can be mapped into the 
square [r, 1 — r] 2 by a translation in the torus. Embeddable components do not 
wrap around the torus. Components which are not embeddable must have a large 
size (at least f2(l/r)). Sometimes several non-embeddable components can coexist 
together. However, there are some non-embeddable components which are so spread 
around the torus that do not allow any room for other non-embeddable ones. Call 
these components solitary. Clearly, by definition we can have at most one solitary 
component. We cannot disprove the existence of a solitary component, since with 
probability 1 — o(l) there exists a giant component of this nature. For components 
which are not solitary, we give asymptotic bounds on the probability of their existence 
according to their size. 

For a fixed integer £ > 1, let Ki be the number of components in G(X\r) of size 
exactly I. For any fixed e > 0, let K' » be the number of components of size exactly 

I which have all their vertices at distance at most er from their leftmost one. Let Kg 
denote the number of components of size > t and which are not solitary. 

Notice that K ' p < Kg < Kg. In the next theorem, we show that asymptotically 

all the weight in the probability that Kg > comes from components which also 
contribute to K' e „ for e arbitrarily small. This implies that at r c , the more common 
components of size > I are cliques of size exactly I, with all their vertices close 
together. 
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Lemma 2.1. Let I > 2 be a fixed integer, and < e < 1/2 also fixed. Assume that 
p = 0(1). Then, 

EK' e/ = 6(1/ log^n) 



Proof. First observe that with probability 1, for each component T which contributes 
to K' p, r has a unique leftmost vertex Xi and the vertex Xj in T at greatest distance 
from Xi is also unique. Hence, we can restrict our attention to this case. 

Fix an arbitrary set of indices J C {1, . . . , n} of size \J\ = I, with two distinguished 
elements i and j. Denote by y = Ufce J ^ k ^ ne se * °^ random points in X with indices 
in J. Let £ be the following event: All vertices in y are at distance at most er from Xi 
and to the right of Xf, vertex Xj is the one in 3^ with greatest distance from Xf, and 
the vertices of y form a component of G{X; r). If P(£) is multiplied by the number 
of possible choices of i, j and the remaining t — 2 elements of J, we get 



In order to bound the probability of £ we need some definitions. Let p = d(Xi,Xj) 
and let S be the set of all points in the torus [0, l) 2 which are at distance at most r 
from some vertex in y. (Notice that p and S depend on the set of random points 3^)- 
We first need bounds of Area (S) in terms of p. Observe that S is contained in the 
circle of radius r + p and center Xi, and then 



Now let i\_ = i, z'r, ij and i& be respectively the indices of the leftmost, rightmost, 
topmost and bottommost vertices in 3? (some of these indices possibly equal). Assume 
w.l.o.g. that the vertical length of 3* (i.e. the vertical distance between Aj T and Xi B ) 
is at least p/y/2. Otherwise, the horizontal length of y has this property and we can 
rotate the descriptions in the argument. The upper halfcircle with center Aj T and 
the lower halfcircle with center Aj B are disjoint and are contained in S. If Aj R is 
at greater vertical distance from Xi T than from X{ B , then consider the rectangle of 
height p/(2 v / 2) and width r - p/{2^/2) with one corner on Xi R and above and to the 
right of Xi R . Otherwise, consider the same rectangle below and to the right of Xi R . 
This rectangle is also contained in S and its interior does not intersect the previously 
described halfcircles. Analogously, we can find another rectangle of height p/(2\/2) 
and width r — p/(2y/2) to the left of Aj L and either above or below Xi L with the same 
properties. Hence, 




(1) 



Area(cS) < vr(r + p) 2 . 



(2) 




(3) 



From ([2]), ([3]) and the fact that p < r/2, we can write 




(4) 
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Now consider the probability P that the n — I vertices not in y lie outside S. Clearly 
P = (1 — Area(5)) n ~ £ . Moreover, by (Jlj) and using the fact that e~ x ~ x < 1 — x < e~ x 
for all x € [0, 1/2], we obtain 

„ „— (l+p/(6r))7rr 2 n 

-(l+5p/(2r))7rr 2 n-(97rr 2 /4) 2 n p e 

(1 - 9vrr 2 /4)^ ' 

and after a few manipulations 

P\ 1+W(2r) e -(9- 2 /4) 2 n p 1 f5l 

Event £ can also be described as follows: There is some nonnegative real p < er 
such that Xj is placed at distance p from X{ and to the right of Xi\ all the remaining 
vertices in y are inside the halfcircle of center X^ and radius p; and the n — £ vertices 
not in y lie outside S. Hence, P(£ ) can be bounded from above (below) by integrating 
with respect to p the probability density function of d(Xi,Xj) times the probability 
that the remaining i — 2 selected vertices lie inside the right halfcircle of center Xi 
and radius p times the upper (lower) bound on P we obtained in ©: 

9(1) 7(5/2) <P(f)< 9(1) 7(1/6), (6) 

where 

m=l *p(- 2 p 2 ) ^r r d P . 



JO 

2 fir , 
- -r' 
n V2 



x 2t-3 n ~Px dx ( 7 ) 



Since £ is fixed, for /3 = 5/2 or /3 = 1/6, 

, log*" 1 
7(/5) = 9 1 ^ 




2£-2 

„ n l log n , 

and the statement follows from (TjQ), ([6]) and ([8]). □ 

Lemma 2.2. 7e£ £ > 2 be a fixed integer. Let e > be also fixed. Assume that 
p = 9(1). Then 

P(K e -K'>0) = O(l/log £ n) 
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Proof. We assume throughout this proof that e < 10~ 18 , and prove the claim for this 
case. The case e > 1CT 18 follows from the fact that (JQ — K' e e ) < (Kg — K' 1Q _ 1S e ). 

Consider all the possible components in G(X;r) which are not solitary. Remove 
from these components the ones of size at most I and diameter at most er, and denote 
by M the number of remaining components. By construction Kg — K' » < M, and 

therefore it is sufficient to prove that P(M > 0) = 0(l/log n). The components 
counted by M are classified into several types according to their size and diameter. 
We deal with each type separately. 

Part 1. Consider all the possible components in G(X; r) which have diameter at most 
er and size between £+1 and logn/37. Call them components of type[H and let A^j] 
denote their number. 

For each k, £ + 1 < k < logn/37, let E\. be the expected number of components 
of type [U and size k. We observe that these components have all of their vertices at 
distance at most er from the leftmost one. Therefore, we can apply the same argument 
we used for bounding EK' e£ in the proof of Lemma I2TT1 Note that (JTJ) , ([6]) and (JZJ) 
are also valid for sizes not fixed but depending on n. Thus we obtain 

E fc <0(l)n(n-l)("~f) 1(1/6), 



k-2, 

where 1(1/6) is defined in ©. We use the fact that < (^) fc_2 and get 



Eh = 0(l)logn (j^)~ 2 f J^n*'* dx (9) 

The expression x 2k ~ s n~ x ^ 6 can be maximized for x E M + by elementary techniques, 
and we deduce that 



2k -3 x 2A '' 



(e/6) log n 

Then we can bound the integral in ([9]) and get 



E k = 0(1) log n 



t logn\ fe 2 / 2k -3 N 2A 



2k -2) \ (e/6) log n 



~ U{l) \2e(k-2)\ogn) * 

/ (2k— 3) 2 \ 

Note that for k < logn/37 the expression ( || (fc_2)iogn ) ^ * s decreasing with k. 
Hence we can write 

Ek = 0[ 1 ), Vfc : ^ + 3< fc< -Uogn. 
\log + n/ 37 

Moreover the bounds Et+\ = 0(1/ log n) and = 0(1/ log £+1 re) are obtained 

from Lemma |2. 11 and hence 



■ log n 

£ 



EMn = V E k = Ol^)+o(—^) + l ^o(—^ r -)=o(^ l 
^ Uog'n/ \\oi +1 n) 37 Vlog^+W \log e 
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and then P(A^J> 0) < EA£j]= 0(l/log £ re). 

Part 2. Consider all the possible components in G(X; r) which have diameter at most 
er and size greater than log re/37. Call them components of typeEl and let denote 
their number. 

We tessellate the torus with square cells of side y = L( er ) _1 J 1 {y — er but also 
y ~ er). We define a box to be a square of side 2y consisting of the union of 4 cells of 
the tessellation. Consider the set of all possible boxes. Note that any component of 
type [2] must be fully contained in some box. 

Let us fix a box b. Let W be the number of vertices which are deployed inside b. 
Clearly W has a binomial distribution with mean EW = (2y) re ~ (2e) 2 logn/vr. By 

log n 

WrEw 



setting 5 = i?-!^ — 1 and applying Chernoff inequality to TV, we have 



1 / e s \ EW (i°g(i+«)-TT7) 
P(W > - log n) = P(W > (1 + 5)EW) < = re 37 . 

Note that 5 ~ j^i — 1 > e 79 , and therefore 

P(W > ^logn) < re -2 " 1 . 

Then taking a union bound over the set of all @(r _1 ) = 0(re/logre) boxes, the prob- 
ability that there is some box with more than ^ log n vertices is (9(l/(n 1 ' 1 log re)). 
Then since each component of type is contained in some box, we have 

P(il^>0)=O(l/(re L1 logre)). 

Part 3. Consider all the possible components in G{X\ r) which are embeddable and 
have diameter at least er. Call them components of type[3j and let A^] denote their 
number. 

We tessellate the torus into square cells of side ar, for some a = a(e) > fixed 
but sufficiently small. Let T be a component of type El Let S = Sr be the set of 
all points in the torus [0, l) 2 which are at distance at most r from some vertex in T. 
Remove from S the vertices of T and the edges (represented by straight segments) 
and denote by S' the outer connected topologic component of the remaining set. By 
construction, S' must contain no vertex in X. 

Now let ii, zr, ij and i& be respectively the indices of the leftmost, rightmost, 
topmost and bottommost vertices in T (some of these indices possibly equal). Assume 
w.l.o.g. that the vertical length of T (i.e. the vertical distance between JQ T and Xi B ) 
is at least er/v2- Otherwise, the horizontal length of T has this property and we can 
rotate the descriptions in the argument. The upper halfcircle with center Xj T and 
the lower halfcircle with center JQ B are disjoint and are contained in S' . If JQ R is 
at greater vertical distance from Xi T than from Xj B , then consider the rectangle of 
height er/(2^2) and width r- er/ (2^2) with one corner on Xj R and above and to the 
right of Xi R . Otherwise, consider the same rectangle below and to the right of Xi . 
This rectangle is also contained in S' and its interior does not intersect the previously 
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described halfcircles. Analogously, we can find another rectangle of height er/(2y/2) 
and width r — er/(2y/2) to the left of X it _ and either above or below X it _ with the same 
properties. Hence, taking into account that e < 10 -18 , we have 



Area(s0 >„ 2+2 (_^) (, 



er \ . e\ 2 



2V2J V 5 



> 1 + - %r\ (10) 



Let 5* be the union of all the cells in the tessellation which are fully contained in S'. 
We loose a bit of area compared to S'. However, if a was chosen small enough, we can 
guarantee that 5* is topologically connected and has area Area(5*) > (1 + e/6)irr 2 . 
This a can be chosen to be the same for all components of typeO 

Hence, we showed that the event (A^gj > 0) implies that some connected union of 
cells S* of area Area(5*) > (1 + e/6)7rr 2 contains no vertices. By removing some cells 
from S* , we can assume that (l + e/6)vrr 2 < Area(S*) < (l + e/6)vrr 2 + aV. Let S* be 
any union of cells with these properties. (Note that there are 0(l/r 2 ) = 0(n/logra) 
many possible choices for S* .) The probability that S* contains no vertices is 



(1 - Area(S*)) n < e - (1+£/6) " ! " = (J^j 



Therefore, we can take the union bound over all the 0(n/logn) possible S* , and 
obtain an upper bound of the probability that there is some component of the type[3j 

*(m> 0) < e (JL.) = e ( . 

'~ \lognJ\nJ \n e / 6 lognJ 

Part 4. Consider all the possible components in G(X;r) which are not embeddable 
but not solitary either. Call them components of type HJ and let A^] denote their 
number. 

We tessellate the torus [0, l) 2 into 0(n/logn) small square cells of side length ar, 
where a > is a sufficiently small positive constant. 

Let r be a component of type[H Let S = Sr be the set of all points in the torus 
[0, l) 2 which are at distance at most r from some vertex in V. Remove from S the 
vertices of T and the edges (represented by straight segments) and denote by S' the 
remaining set. By construction, S' must contain no vertex in X. 

Suppose there is a horizontal or a vertical band of width 2r in [0, l) 2 which does 
not intersect the component V (assume w.l.o.g. that it is the topmost horizontal band 
consisting of all points with the y-coordinate in [1 — 2r, 1)). Let us divide the torus 
into vertical bands of width 2r. All of them must contain at least one vertex of T, 
since otherwise T would be embeddable. Select any 9 consecutive vertical bands and 
pick one vertex of V with maximal y-coordinate in each one. For each one of these 9 
vertices, we select the left upper quartercircle centered at the vertex if the vertex is 
closer to the right side of the band or the right upper quartercircle otherwise. These 
nine quartercircles we chose are disjoint and must contain no vertices by construction. 
Moreover, they belong to the same connected component of the set S' , which we denote 
by S", and which has an area of Area(5 // ) > (9/4)7rr 2 . Let S* be the union of all the 
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cells in the tessellation of the torus which are completely contained in S" . We lose a 
bit of area compared to S" . However, as usual, by choosing a small enough we can 
guarantee that S* is connected and it has an area of Area(5*) > (ll/5)7rr 2 . Note that 
this a can be chosen to be the same for all components T of this kind. 

Suppose otherwise that all horizontal and vertical bands of width 2r in [0, l) 2 
contain at least one vertex of T. Since T is not solitary it must be possible that it 
coexists with some other non-embeddable component T' . Then all vertical bands or 
all horizontal bands of width 2r must also contain some vertex of V (assume w.l.o.g. 
the vertical bands do). Let us divide the torus into vertical bands of width 2r. We can 
find a simple path II with vertices in r' which passes through 11 consecutive bands. 
For each one of the 9 internal bands, pick the uppermost vertex of T in the band 
below II (in the torus sense). As before each one of these vertices contributes with 
a disjoint quartercircle which must be empty of vertices, and by the same argument 
we obtain a connected union of cells of the tessellation, which we denote by S* , with 
Area(5*) > (ll/5)7rr 2 and containing no vertices. 

Hence, we showed that the event (A^j > 0) implies that some connected union 
of cells S* with Area(5*) > (ll/5)7rr 2 contains no vertices. By repeating the same 
argument we used for components of type[3]but replacing (1 + e/6)irr 2 by (ll/5)-7rr 2 , 
we get 

p (^>») = e(^5^). 

□ 

Lemma 2.3. Let £ > 2 be a fixed integer. Let < e < 1/2 be fixed. Assume that 
H = 0(1). Then 

E[<,] 2 = 0(l/log 2 ^ 2 n) 

Proof. As in the proof of Lemma 12.11 we assume that each component V which con- 
tributes to K' ei has a unique leftmost vertex Xi, and the vertex Xj in T at greatest 
distance from Xi is also unique. In fact, this happens with probability 1. 

Choose any two disjoint subsets of {1, . . . ,n} of size t each, namely J\ and J2, 
with four distinguished elements £ J\ and 12, 32 £ 'h- For k £ {1,2}, denote by 
y k = |J ig j fc X[ the set of random points in X with indices in J^. Let £ be the event 
that the following conditions hold for both k = 1 and k = 2: All vertices in y k are 
at distance at most er from X,- lk and to the right of Xi k ; vertex Xj k is the one in 
with greatest distance from Xi k ; and the vertices of 34 form a component of G{X; r). 
If P(£) is multiplied by the number of possible choices of ifc, jt and the remaining 
vertices of we get 

E[K> ie ]2 = 0(n 2 ")P(£). (11) 

In order to bound the probability of £ we need some definitions. For each k G 
{1,2}, let pk = d(X-i k , Xj k ) and let Sk be the set of all the points in the torus [0, l) 2 
which are at distance at most r from some vertex in y^. (Obviously and depend 
on the set of random points 3V) Also define S = Si U 1S2. 
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Let T be the event that d(Xi 1 , Xi 2 ) > 3r. This holds with probability 1 — 0(r 2 ). 
In order to bound P(£ | J 7 ), we apply a similar approach to the one in the proof of 
Lemma [2. II In fact, observe that if T holds then S\ HS2 = 0. Therefore in view of ([!]) 
we can write 

vrr 2 (2 + ( Pl + p 2 )/(6r)) < Area(S) < ^V, (12) 
and using the same elementary techniques that gave us ([5]) we get 

, m ,. of fll\ 2 HPi+P2)/<S>r) 1 

(1 - Area(S))»-- < (f) (1 _ 18)rrW r (13) 

Now observe that £ can also be described as follows: For each k £ {1,2} there is some 
nonnegative real pt < er such that Xj. is placed at distance pk from X{ k and to the 
right of Xi k ; all the remaining vertices in 34 are inside the halfcircle of center Xi k 
and radius p^', and the n — I vertices not in 34 he outside Sk- In fact, rather than 
this last condition we only require for our bound that all vertices in X\ (34 U3^) are 
placed outside S. Clearly, this has probability (1 — Area(5)) n_2 ^. Then, from (| 13j) and 
following an analogous argument to the one that leads to ([6]), we obtain the bound 

P(S I T) < 9(1) / / vrpi ( -p 2 ) vrp 2 ( -p 2 ) d Pl dp 2 
Jo Jo 



2 H1J "^V2 F2 / n 2+(p 1+ p 2 )/(6r) 

2 



= 9(1)1(1/6) 

where 1(1/6) is defined in ©. Thus from (jHJ) we conclude 

p(£ a < e(i) p{F) /(i/6) 2 = o ( n2qog V 2 J ■ d4) 

Ctherwise, suppose that T does not hold (i.e. d(Xi 1 , Xi 2 ) < 3r). Cbserve that S 
implies that d{Xi 1 ,Xi 2 ) > r, since and Xi 2 must belong to different components. 
Hence the circles with centers on X^ and Xi 2 and radius r have an intersection of 
area less than (7r/2)r 2 . These two circles are contained in S and then we can write 
Area(5) > (3/2)7rr 2 . Note that £ implies that all vertices in X \ (34 U 34) are placed 
outside S and that for each k £ {1, 2} all the vertices in 34 \ {Xi k } are at distance at 
most er and to the right of Xi k . This gives us the following rough bound 

Multiplying this by P(T) = 0(r 2 ) = O (log n/n) we obtain 

^^) = 0[ l ^f], (15) 



n 



which is negligible compared to (114D . The statement follows from (jlll) , (I14j) and fjl5[) . 

□ 
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The main result of this section now follows easily. 

Theorem 2.2. Let I > 2 be a fixed integer. Let < e < 1/2 be fixed. Assume that 
[i = 0(1). Then 



Pr 



K, > 



Pr > 0] ~ Pr [K' e/ > 0] = 6 (j^^) 



Proof. From Corollary 1.12 in [2|, we have 

EK' e>e - < P«* > 0) < EK> e/ , 

and therefore by Lemmata 12.11 and 12.31 we obtain 

P(<,>0) = 6(l/log^ 1 n). 
Combining this and Lemma 12.21 yields the statement. □ 

3 Dynamic Properties 

We define the dynamic model as follows. Given a positive real s = s(n) and a 
positive integer m = m(n), we consider the following random process (Xt)t& = 
(X t (n, s,m)) t <=z: At time step t = 0, n agents are scattered independently and u.a.r. 
over the torus [0, l) 2 , as in the static model. Moreover each agent chooses u.a.r. an 
angle a G [0, 2ir), and moves in the direction of a, travelling distance s at each time 
step. These directions are changed every m steps for all agents. More formally, for 
each agent i and for each interval [t, t + m] with i 6 Z divisible by m, an angle in 
[0, 2ir) is chosen independently and u.a.r., and this angle determines the direction 
of i between time steps t and t + m. Note that we are also considering negative 
steps, which is interpreted as if the agents were already moving around the torus ever 
before step t = 0. We extend the notation from the static model, and denote by 
Xi t t = {xi : t,yi,t) the position of each agent i at time t. Also let X t = (JiLi -^*>< ^ e ^ ne 
set of positions of the agents at time t. Furthermore, given a positive r = r(n) G R 
such that r = o(l), a random graph process can be derived from (Xt)tel,- F° r an y 
t G Z, the vertex set of G{Xt;r) is X t , and we join by an edge all pairs of vertices in 
Xt which are at Euclidean distance (in the torus) at most r. We derive asymptotic 
results on [G(Xt;r)) t 2 as n — > oo. 

We use the following lemma proven in |llj . 

Lemma 3.1. At any fixed step t e Z, the vertices are distributed over the torus 
[0, l) 2 independently and u.a.r. Consequently for any t £ Z, G(Xt\r) has the same 
distribution as G(X;r). 

In the remaining of the section, we focus our attention around the threshold of 
connectivity and we assume that fi = 0(1), or equivalently 



/logndbO(l) 



nn 
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In order to prove the main statement of this section, we first require some technical 
results which involve only two arbitrary consecutive steps t and t+1 of {Xt)t^z- In this 
context t is considered to be an arbitrary fixed integer, and it is often omitted from 
notation whenever it is understood. Thus for each i G {1, . . . , n}, the random positions 
Xi t and -Xj,t+i of agent i at times t and t + 1 are simply denoted by X{ = (xi,t/i) and 
X[ = (x'^y'j). Let also X = X t and X' = X t+ i. Note that the random points Xj and X[ 
are not independent. In fact if 2nzi (zi G [0, 1)) is the angle in which the agent i moves 
between times t and t + 1, then x\ = Xi + s cos(2-KZi) and y[ = Ui + s sin(27TZj) (where 
all the sums involving coordinates are taken mod 1). This motivates an alternative 
description of the model at times t and t + 1 in terms of a three-dimensional placement 
of the agents, in which the third dimension is interpreted as a normalized angle. For 
each i G {l,...,n}, define the random point Xi = (xi,yi,Zi) G [0, l) 3 , and also let 
X = UiLi^*- Observe that by Lemma 13.11 all the random points Xj are chosen 
independently and u.a.r. from the 3-torus [0, l) 3 , and also that X encodes all the 
information of the model at times t and In fact, if we map [0, l) 3 onto [0, l) 2 

by the following surjections 

7Ti : [0,1) 3 ^ [0,1) 2 vr 2 : [0, l) 3 - [0, l) 2 

(x, y, z) i — > (x, y) (x, y,z) \-* (x + s cos(2nz),y + s sin(27T2;)), 

we can easily recover the positions of agent i at times t and t + 1 from Xj and 
write Xi = 7Ti(Xj) and X[ = ir^Xj). Notice moreover that, for any measurable set 
•A C [0, l) 2 , the events Xi G A and X[ G A are respectively equivalent to the events 
Xi G 7rj~ 1 (^4) and X{ G 7r^ (A) in this new setting. Furthermore, we have 

Area(^) = Vo\fc\A)) = Vcifc 1 (•*)) , (16) 

since 

Vol^ 1 ^)) = Vo\(A x [0, 1)) = Area(^), 
and also by putting A z = A — (scos(2ttz), ssin(27rz)) 

Vo\(tt^ 1 (A))= [ ([ dxdy)dz = Area(A). 
J [0,1) v Ja z ' 

Naturally (|16() is compatible with the fact that, in view of Lemma l3.1( for any mea- 
surable sets A C [0, l) 2 and B C [0, l) 3 , 

P(Xi eA) = Area(^4), P(X^ G A) = Area(^), P(X G B) = Vol(B). 

Now we define some sets which will repeatedly appear in this section. For each 
i G {1, . . . , n}, consider the sets 

Ki = {X G [0, l) 2 : d(X,Xi) < r} and % = {X G [0, l) 2 : d{X,X[) < r}, 

and also let Hi = 7rj~ 1 (7?.j) and 1Z[ = 7r^~ 1 (7^ / i ) be their counterparts in [0, l) 3 . Note 
that for i,j G {l,...,n}, we have that Xj G IZj iff d(Xi,Xj) < r, and similarly 
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that Xi G IZ'j iff d(X 4 ',Xj) < r (each of these events occurring with probability 
exactly Vol(7£j) = Vol^'j) = irr 2 ). Also observe that Xj is isolated in G(X;r) iff 
(X\{Xi})nKi = 0, and that analogously X[ is isolated in G(X'; r) iff (,£\{Xj})nft- = 
0. 

We need the following 

Lemma 3.2. Assume \i = 0(1). There exists a constant e > swc/i i/iai the following 
statements are true (for large enough n): For any i,j G {1, . . . , n} (possibly i = j), 

1. ifd(Xi,Xj) > r then Vol(7^ n 7^) < § r 2 . 

2. i/a < r/7 anddiXi^Xj) >r-2s then Vol((7^ U 7^) n U 7^)) < (l-e)vrr 2 . 
5. i/s > r/7 and s = O(r) f/ien Vol(7^ n7^) < (1 - e)7rr 2 . 

^. ifs = u(r) then Vol (7^n7^) = 0(r 3 ^±±) = o(r 2 ). 
Proof. 

Statement 1. Assume w.l.o.g. that the segment XjXj is vertical and that Xj is above 
X,-. Then, let 5 C [0, l) 2 be the upper halfcircle with center Xj and radius r, and 
5 = ^(S) = S x [0, 1) C [0, l) 3 . Clearly, Vol(S) = vrr 2 /2, ScRj and 5nK i = 0, 
and the statement follows. 

Statement 2. The distance between X 4 ' and Xj is greater than 3r/7, since d(X|, X 1 -) > 
d(Xi,Xj) — 2s > r — 4s. Let Sj (respectively Sj) be the set of points in [0, l) 2 at 
distance at most 8r/7 from X[ (respectively Xj) Note that <Sj and Sj are two circles 
of radius 8r/7 with centers at distance greater than 3r/7. Then straightforward 
computations show that Area(Sj r\Sj) is at most (1 — e)7rr 2 for some e > 0. We define 
Si = 7r^~ 1 (»Si ) and <Sj = 7r^ 1 (6'j). Clearly, Si D 7£j U T^j and Sj D 7^.j U 7£j. Hence, 

Vol((£i U ft-) n (Wj- U Wj-)) < Vol(5i n = Area(Sj n 5,-) < (1 - e)nr 2 . 

Statement 3. Let k G {1, . . . , n} be different from i and j. Observe that \/o\(1Zi\1Z'j) is 
the probability that d(Xj, X&) < r but d(X'- , X£) > r. Suppose that d(Xj, X&) < r but 

also d(Xj,X fc ) > 13r/14. (This happens with probability at least (1 - 13 2 /14 2 )vrr 2 .) 

> 

Let a be the angle of X'-X^ with respect to the horizontal axis. Recall that agent k 
moves between time steps t and t + 1 towards a direction 2irzk, where zt is the third 
coordinate of X&. If 2-7^ G [a— 7r/3, a+ir/S], then the agent increases its distance with 
respect to X'- by at least s/2 > r/U and thus d{X' j ,X' k ) > r/U + 13r/14 = r. This 

range of directions has probability 1/3. Summarizing, we proved that Vol(7£; \Tl'j) > 
(1 - 13 2 /14 2 )7rr 2 /3, and the statement follows. 

Statement 4. Given k G {1, . . . ,n} different from i and j, observe that Vol(7£j n TZj) 
is the probability that d(Xfc,Xj) < r and also d{X' k ,X'-) < r. Suppose first that 
s < 1/2. We claim that the probability that d(X' k ,X'j) < r conditional upon any 
fixed outcome of X& is at most (2 + e)r/s for some e > 0, no matter which particular 
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point Xk is chosen. In fact, assume Xk ^ X'- (the case Xk = X'- is trivial) and let a 

> 

be the angle of X^X'j with respect to the horizontal axis. If agent k moves between 
steps t and t + 1 towards a direction 2irzk not in [a — arcsin(r/s), a + arcsin(r/s)] 
then d(X k ,X'j) > r. Hence, Vol (7^ n 7^) is at most P(d(X k ,Xi) < r) = irr 2 times 
(2 + e)r/s, which satisfies the statement. 

The case s > 1/2 is a bit more delicate, since agent k may loop many times 
around the torus while moving between steps t and t + 1. In fact, as we move along 
the circumference of radius s centered on Xk we cross the axes of the torus 0(1 + s) 
times. This gives the extra factor (1 + s) in the statement, which is negligible when 
s = o(l) but grows large when s = u;(l). □ 

For each % G {1, . . . , n}, we define Qi = TZ' t \ IZi and Q[ =1Zi\lZ' i . Given any two 
agents i and j, observe that Xj G Q'j iff Xj G Q- iff d(Xi,Xj) < r and d(X' i ,X' j ) > r 
(i.e. the agents are joined by an edge at time t but not at time t + 1). This holds with 
probability Vol(Qj) = Vol(Q^), which does not depend of the particular agents and of 
t and will be denoted by q hereinafter. The value of this parameter depends on the 
asymptotic relation between r and s and is given in the following 

Lemma 3.3. The probability that two different agents i,j G {1, . . . ,n} are at distance 
at most r at time t but greater than r at time t+1 is q < nr 2 , which also satisfies 

^sr ifs = o(r), 
G(r 2 ) ifs = Q(r), 
7rr 2 if s = io(r). 

Proof. The first bound on q is immediate from the definition of q and the fact that 
Vol(Tvj) = irr 2 . In order to obtain the second statement, we consider separate cases. 

Case 1 (s < er, for some fixed but small enough e > 0). In order to compute the 
probability that Xj G Q[, we express Xj = (xj,yj,Zj) in new coordinates (p,6,z), 
where p = d(Xj,X^), is the angle between the horizontal axis and X{Xj, and 
z = Zj. Then we integrate an element of volume over the region in terms of these 
coordinates. Let us also call £ = d(Xj,Xi), so that (£, 6,z) are the usual cylindrical 
coordinates. From the law of cosines, we can stablish the relation between p and £ 
and write 



p = _|_ s 2 _ 2£scos6> and £ = J p 2 - s 2 sin 2 6 + s cos 9. (17) 

Now observe that the minimum value that p can take is r — s, since Xj must lie 
outside the circle of radius r — s and center X[ . Otherwise by the triangular inequality 
d(X' i , Xj) < r and the agents i and j would share an edge at step t + On the other 
hand, Xj must lie inside the circle of radius r centered on JQ, and therefore (by setting 
£ = r in (|17j) ) the maximum value that p can achieve is 

p = \J r 2 + s 2 — 2rs cos 6 
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Moreover, let a be the angle determined from the range of all possible values of 2-kz 
(i.e. possible directions for agent j to move). Again by the law of cosines, 



a = 2 arccos 



r 2 — s 2 — p 2 
2sp 



Finally from (|17p and the change of variables formula, it is straighforward to determine 
the element of volume in coordinates (p, 9, z): 

dxdydz = £ d^dOdz = - — — — - dpd9dz, 
c, — s cos 9 

where, using the fact that r — 2s < £ < r, we can write 



£ — s cos 

In view of all the above, we deduce 

q = / dxdydz 
JQ! 

2ty / t y/r 2 +s 2 —2rs cos 6 p 

— - ^ - dpde 

2ir £ — s cos 

/ ti \ \ f2ir r\/r' 2 -s 2 sin 2 <f>-s cos tf> i / 2 _ „2 _ 2\ 

1±0 (;))i L k^—hr-)^ 
-( 1±0 (i}) 2 l"U-" Bme - er2 

9 2 n\ r COS 6 — S \ 

+ (r +s — 2rs cos arccos w. 

V r 2 + s 2 — 2rs cos ^ / 

Now by looking at the Taylor series with respect to s/r of the expression inside the 
integral divided by r 2 , we get 

s \ \ n 2 ( 29 cos 9 s ^ ( ( s 



1±0 UU r {- — r + {lr) )) dl>= ( 1±0 {-r))^ r - (18) 

Case 2 (er < s < r/7). Recall that IZi is the circle of radius r and center JQ. 
Take the chord in 7£j which is perpendicular to the segment XiX- and at distance 
r from This chord divides TZi into two regions. One of them (call it S) has the 
property that all the points inside are at distance at least r from X[ and moreover 
Area (S) > e\/2e — e 2 r 2 . Suppose that Xj G S (i.e. the agent j is in S at time t), 
which happens with probability at least e\/2e — e 2 r 2 . Let us now consider the circle 
centered on X[ and passing through Xj. We observe that d(X'j,X' i ) > d(Xj,X' i ) 
with probability at least 1/2, since it is sufficient that the direction 2irzj in which 
agent j moves lies in the outer side of the tangent of that circle at Xj. Therefore, 
the probability that d(Xj,Xi) < r and d(X'j,X' i ) > r (i.e. Xj S is at least 
iev / 27^V 2 . 
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Case 3 (s > r/7). We can write 

q = Vol(Qj) = Vol(^i \ ft-) = Vol(^) - Vol(^j n ft-), 
and the result follows from the statements 3-4 in Lemma 13.21 □ 
We also need the following 

Lemma 3.4. Consider a setting with n balls and k + 1 disjoint bins Uq, . . . ,Uk, where 
each ball is either placed into at most one of the bins or possibly remains outside all 
of them. Suppose that each ball is assigned to bin Ui with probability pi = Pi(n), 
independently from the choices of the other balls. Moreover suppose that for all i 
(0 < i < k) we have pi = o(l), where asymptotics are with respect to n — > oo and k 
is assumed to be fixed. Then, the probability P that Uq contains no balls but each Ui 
has at least one for all 1 < i < k is 



p~(i- P0 rii(i-e- n n- 



1=1 



Proof. By reordering the labels of the bins except for Uo, assume that npi = o(l) if 
1 < i < r and npi = 0(1) if r + 1 < i < k. For any / = o(l), we define 



Pf= E H) EUai (i-%r)" 



ai,...,a r 

where the summation indices a\ , . . . , a r run from to 1 . We can write 



Pf 



ai,...,a r mo,...,m r >0 
moH \-m r =n 



e H> 5 " e L'jnr 



mo, . . .,m r J 11 \ l - f 



with the convention 0° = 1. A changing of the order of summation converts this 
expression to 



E >o L,.".,Jn(3) EMF-n- 

,...,77l r >0 

moH Vvn r =n 



.mn. , 

mo,...,m r >0 i=l ai,...,a r i=l 



r 

m,i 
i 



All terms in this sum with some nii = (1 < i < r) cancel, since each of these terms 
is equal but has opposite sign to the one obtained by switching the value of aj. So, 
only the terms with all mi, ■ ■ ■ , m r > 1 remain, and among these we can remove the 
ones with some en = 0. Hence, 



]?>o ( ir U,---,mjn(i Pi 

mi...,m r >l 
mo-\ \-m r =n 
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Since npi = o(l), the main asymptotic weight in this sum corresponds to the term 
mo = n — r and mi, . . . , m r = 1, so 



P f ~ TTZTF II« ~ Eh* ~ II( 1 - e " nPl )- (19) 
^ J ' i=l 1=1 i=l 

By an inclusion-exclusion argument, the probability in the statement can be written 

as 

(k \ n 
i-po-E'Hw) ' 
1=1 / 

where the summation indices oi, . . . , at run from to 1. Then, if we define 



E[=i (H Pi 

a r +i,...,a k = y-L)— -ix- — r 



ai,...,a 



E (-i) EUai i-- 



we can write 



P =(!-«»)" E (-1)^— fl-^^) P ar+ll ...,a 



a r +i,...,aj. 



1 -Po 

= (l-Po) n E (-ll^^pf-fl + oflDEwi]^ (20) 

a r +i,...,a fc \ i=r+l / 

Note that for each a r+ i, . . . , S {0,1}, in view of (fT9j) and setting f = Po + 

YH=r+l a iPi> We have 

r 

P ar+1 ,...,a, ~ II( 1 " e " nPl )- ( 21 ) 
i=l 

The fact that npi = $7(1) for r + 1 < i < m prevents the leading term of the sum 
in (|20p from cancelling out. Thus, from (|20p and (|2ip . we obtain 



p~(i-po) n n( i - e_npi )- ° 

i=l 

We are now in good position to study the changes experienced by the isolated 
vertices between two consecutive steps t and t+1. Extending the notation in Section[21 
we denote by Kit the number of isolated vertices of G{Xt,r). Also, for any two 
consecutive steps t and t + 1, we define the following random variables: Bt is the 
number of agents i such that Xi is not isolated in G(Xt;r) but X[ is isolated in 
G(Xt+i; r); Dt is the number of agents i such that Xi is isolated in G(Xt; r) but X[ is 
not isolated in G(Xt+i;r); St is the number of agents i such that Xi and X[ are both 
isolated in G(X t ;r) and G(X t+ i;r). For simplicity, we often denote them just by B, 
D and S whenever t and t + 1 are understood. Note that B and D have the same 
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distribution, since any creation of an isolated vertex corresponds to a destruction of 
an isolated vertex in the time-reversed process and vice versa. 

To state the following result, we need one more definition: Given a collection of 
events £i(n), . . . , £fc(n) and of random variables Wi(n), . . . , Wi(n) taking values in N, 
with k and I fixed, we say that they are mutually asymptotically independent if for 
any k', I', i\, . . . , ik',ji, ■ ■ ■ , jy, • • • > W V e N such that k' < k, I' < I, 1 < i\ < ■ ■ ■ < 
iy < k, 1 < j\ < ■ ■ ■ < ji> < I we have that 

(k' V \ k' V 

/\£ ia Af\ (W Jb = w b ) ~ fj P(£ la ) J] *{W jb = w b ). (22) 
a=l 6=1 / a=l 6=1 

Proposition 3.1. Assume \x = Then for any two consecutive steps, 

EB = ED ~ n{l - e~ qn ) and ES ~ /ie" 9 ". 
Moreover we have that 

1. If s = o(l/rn), then P(B > 0) ~ E5; P(L> > 0) ~ ED; S is asymptotically 
Poisson; and {B > 0), (D > 0) and S 1 are asymptotically mutually independent. 

2. If s = @(l/rn), t/ien D and S are asymptotically mutually independent 
Poisson. 

3. If s = u(l/rn), then B and D are asymptotically Poisson; P(S > 0) ~ E5; and 
B, D and (S > 0) are asymptotically mutually independent. 

Proof. The central ingredient in the proof is the computation of the joint factorial 
moments E ([B]^ [D]i 2 [S]e 3 ) of these variables. In particular, we find the asymptotic 
values of E (B), E (D) and E (S). Moreover, in the case s = 0(l/(rn)) , we show that 
for any fixed naturals £±, £2 and £3 we have 

E([B]eAD]i 2 [Sh)-^(B^E(DfE(S) e \ (23) 

Then, the result follows from Theorem 1.23 in [2]. The other cases are more delicate 
since ([23]) does not always hold for extreme values of s, and we obtain a weaker result. 
In the case s = o(l/(m)), we compute the moments for any natural £3 but only for 
G {0,1,2} and obtain 

E([B] tl [D} h [S} i3 ) ~ (E(B)*(E (E(S))* if < 2, 
E([B] il [D] i2 [S] l3 ) = o(E(B[D} l2 [S) h )), 

E([B] il [D] i2 [S}e 3 ) = o(E([B] £l D[S} i3 )). (24) 

From this and by using upper and lower bounds given in [2], Section 1.4, applied to 
several variables, we deduce that (B > 0), (D > 0) and S satisfy (|22|) and also 

E ( 9^ 

Pr [B > 0] ~ E (B), Pr [D > 0] ~ Pr [D] and Pr [S = k] ~ e~ E(5) 1 ; VA; G N. 
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Similarly, in the case s = co(l/(rn)), we compute the moments for any naturals £\ 
and £2 but only for £3 G {0, 1, 2} and obtain 

B([B] ei [D] e2 [S] i3 ) ~ (E(B)^(E(D)Y^E(S)Y\ if i 3 < 2, 
E([B] ei [D] e2 [SU 3 ) = o(E([B] h [D] h S)) (25) 

From this and by using once more upper and lower bounds given in Section 1.4 of [2], 
we conclude that B, D and (S > 0) satisfy (f22l) and also 

Pr [B = k] ~ e" E ( g ) (E( 5 )fc VfceN, 

re! 

Pr [£> = fe] ~ e- E ( g ) (E ( ,^ ))fc VreGN and Pr [S > 0] ~ E (S). 

re! 

First, we define for each i G {1, . . . ,n} and as the indicator functions of 

the following events respectively: X$ is not isolated in G(Xt;r) but X[ is isolated in 
G(X t +i;r); Xi is isolated in G(Xt;r) but X[ is not isolated in G{X t+ \\r); Xi and X[ 
are both isolated in G(Xt\r) and G(Af_|_i;r). This allows us to write 

n n n 

B = Y,Bi, D = Y J D l , S = Y,Si- 

i=l i=l i=l 

Note that Bi = 1 iff all points in X \ {Xi} are outside 1Z' { but at least one is inside 
Q[] also Di = 1 iff all points in <Y \ {Xj} are outside IZi but at least one is inside Qf, 
and finally Si = 1 iff all points in X \ {Xi} are outside Hi U % = K i uQ, i = % U 

Now given any fixed naturals £\ , £2 , ^3 with sum £ = l\ + £2 + ^3 , we choose an 
ordered tuple J of I different vertices G {1, . . . , n}, and define 

£=/\(B ta = l)A /\ (A b = l)A /\ (S ic = l). (26) 

a=l fe=^i+l c=h+h+l 

Observe that P(£) does not depend on the particular tuple J, and multiplying it by 
[n]i (i.e. the number of ordered choices of J) we get 

E([B] il [D] l2 [S} i3 ) = [n} i P(£) (27) 

By relabelling the vertices in J we assume hereinafter that J = (1, . . . ,£), and we call 
y = Ui = i{Xi}. Moreover, we define the set 

£=|J£ju(J£iU U (^ U ^) 
i=i i=ii+i i=e 1 +e 2 +i 

and the collection of sets 

Q = {Q'x, ■■■,Q[ 1 , Qti+i, • • • , Qii+e 2 }: 

which play an important role in the computation of P(£)- It is useful to call Q* = Q'i 
for I < i < h, Q* = Qi for £ x + 1 < i < £ x + £ 2 , so that we can write Q = 
{Qi, • • • , Q*i 1+i2 }- 
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Case 1 (s = 0(l/(rn)) ). We say that a vertex % £ J is restricted if there is some 
other j £ J with j > i such that d(Xi,Xj) < 2r + 4s. Let J 7 be the event that there 
are no restricted vertices in J, i.e. d(Xi,Xj) > 2r + 4s for all i, j £ J (i ^ j). This has 
probability 1 — 0(r 2 ). We first suppose that J 7 holds and compute the probability of 
£ conditional upon that. We observe that J- implies that for any i, j £ J (i ^ j) we 
must have TZiDKj = 0, K^nUj = and UiHUj = 0. Then Vol (ft) = iirr 2 +£ 3 q, and 

the sets in Q are pairwise disjoint and also disjoint from ft. Moreover observe that, 
conditional upon J 7 , £ is equivalent to the event that all points in X \ y lie outside 
ft, but at least one belongs to each Q* £ Q. From all the above, the probability of 
£ can be easily obtained by a Balls and Bins argument. In fact, from Lemmata 13.31 
and 13.41 we conclude 

P(£ AT) = (l-0(r 2 ))P(£ | T) 

~ (1 - £vrr 2 - £ 3 q) n (l - e ~ qn ) h+ ^ 

(£.y (1 - e -m)li+ta e -ia9n m ( 28 ) 

We claim that this is the main contribution to P(£). In fact if T does not hold (i.e. 
some of the points in y are at distance at most 2r + 4s), then P(£ | T) is larger 
than the expression in ([25]) . but this is balanced out by the fact that P(.F) is small. 
Before proving this claim, define Ti to be the event that d(Xi,Xj) > r — 2s for all 
i,j £ J (i ^ j). Notice that £ implies H, since otherwise, for some i, j £ J, X; L and 
Xj would be joined by an edge in G{Xt, r) and also X[ and X'- in G{Xt+\\ r), which is 
not compatible with £. Therefore we only need to see that P(£ f\J-) = P(J- AH)P(£ | 
J 7 AH.) is negligible compared to ([28]) . 

Suppose then that Ti. holds and also that p > of the vertices in J are restricted 
(i.e. T does not hold). This happens with probability 0(r 2p ). In this case, we deduce 
that Vo\(TZ) > {I — p)irr 2 + e7rr 2 , since each unrestricted vertex in J contributes at 
least 7rr 2 to Vol (TZ) and the first restricted one gives by Lemma 13.21 (2) the term eirr 2 . 
Moreover, £ implies that all points in X \ y lie outside of TZ, which has probability 
(l — Vo\(TZ)) n ^ = 0(l/n^~ p+e ). Summarizing, the weight in P(£ A J-) coming from 
situations with p restricted vertices is 0(r 2p /n i ~ p+e ) = 0(log p n/n e+e ), and is thus 
negligible compared to ([25]) . Hence P(£) ~ ~P(£ A J 7 ), and the required condition on 
the moments announced in (1230 follows from (j2T[) and (1280 . 

Case 2 (s = o(l/(rn)) ). Defining J 7 and Ti as in the case s = 0(l/(rn)) and by an 
analogous argument we obtain 

P{£ AJ 7 )- (^{l - e -i*)*i+t2 e -t3*n „ ^ ( qn yi+h (29) 

However, the analysis of the case that J 7 does not hold is slightly more delicate here. 
Indeed, there is an additional o(l) factor in (I29p . namely (qnY 1+e2 , which forces us 
to get tighter bounds on P(£ A T A Ti) than the ones obtained before. Unlike in the 
case s = 0(1/ (to)), here we need to consider also the role of Q when J 7 does not 
hold, and special care must be taken with several new situations which do not occur 



21 



otherwise. For instance, since the elements of Q are not necessarily disjoint, then for 
Q*, Q*j G Q the condition that both contain some element of X can be satisfied by 
having just a single point in Q* n Q* fl A'. Moreover, if > 2 and 1 < i < j < l\ (or 
£2 > 2 and £\ + l < i,j < £1 + £2), the previous condition is also satisfied if Xj G Q*, 
which is equivalent to X{ G Q*. If the latter situation occurs, we say that i and j 
collaborate. 

We now distinguish two cases whether two vertices both belonging to J collaborate 
or a vertex from outside J causes T not to hold. First, we bound the weight in P(£AJ-) 
due to situations in which there are no pairs of elements in J which collaborate. We 
need some definitions. Let J\ = {1, . . . ,£\ + £2} and 3^1 = Uiil^ 2 {^?}; an d consider 
the class & of partitions of J\. Namely, a partition of J\ is a collection of blocks 
(nonempty subsets of Ji) which are disjoint and have union J\. The size of a partition 
is the number of blocks, and for each block we call leader to the maximal element in the 
block. Given a partition P = {Ai, . . . , A^} G 2? and also i%, . . . , i}~ G {1, . . . , n}\ J, let 
£p,h,...,i k be the following event: For each block Aj of P, we have G C\ieAj QX an< ^ 
moreover all the points in X\(yu{ii, . . . , ik}) lie outside of 7Z. We wish to bound the 
probability of Sp^,...^ AJ^ATi. Notice that if £p,i Xt ... t i k holds, then all the £\ + £2 — k 
non-leader elements in J\ must be restricted, and possibly some other p' vertices in J 
are restricted too. Moreover, T does not hold iff this p' satisfies < £\+£2~ k+p' < £. 
Given any p' with that property, suppose that the number of restricted vertices in J 
which are either in J\J\ or are leaders of some block is exactly p' . We condition upon 
this and also upon TL, which has probability r 2p . Then for each block Aj with leader 
lj, event £p t i 1: ... t i k requires that G Q*. and for all i G Aj (i 7^ lj) Xi G {Qi 3 U Q^.). 
In addition, since the number of restricted vertices in J is £\ +£2 — k +p' > 0, arguing 
as in the case s = 0(l/(rn)) we have \Zo\(1Z) > (£3 + k — p')7rr 2 + eirr 2 . Then the 
contribution to P(£p t i l! ...,i k A T A H) for this particular p' is 

0(^V( 2 ,)^-*(1 - \,o\{U)Y- e - k = O (on)^ e \ 

so for some < e' < e we can write 

Finally observe that if there are no pairs of elements in J which collaborate, then 
£ A J 7 implies that £p ! i 1 ,... t i k A T A TL holds for some P G 8? of size k and some 
ii, . . . , ifc G {1, . . . , n} \ J, and therefore has probability 

O H O (^) fen)*"". = O (^) ( q n)^, (30) 

negligible compared to (|29p . In particular, if £i,£2 < 2, then no pair of elements in J 
collaborates and then P(£) ~ P(£ A J 7 ). Hence, the first line of (|24l) follows from (|27p 
and ([291). 
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We now extend the approach above to deal with situations in which some pair of 
elements in J collaborate. Unfortunately, their contribution to P(£ A T A TL) may 
be larger than (|29|) if s tends to fast. Hence we restrict £\ and £2 to be at most 2 
and prove only (f2l"|) . If £\ = 2 let £1 be the following event: Xi G Q' 2 , 7£ contains no 
points in \ and also, for each natural (3 < i < 2 + .£ 2 ), Qi contains some point in 
X \ y. Similarly if £ 2 = 2 let £ 2 be the following event: G ^ contains 

no points in X \ y and also, for each natural i (1 < i < £\), Q[ contains some point 
in X \ y. Finally if t\ = £2 = 2 let £ i j2 be the following event: X\ G Q 2 and also 
^3 G Q4. In order to compute P(£i AW), we can repeat the same argument above, 
but imposing that X\ G Q' 2 and ignoring other conditions on Q[ and Q' 2 . We obtain 
for some e' > 

P(£i A W) = O ff(«n) ft = O (an)^, (31) 

and similarly 

P(£ 2 A W) = O (-^) (qn) h+1 and P(£ li2 A H) = O (-^ \{qnf. (32) 

Observe that if some vertices in J collaborate, then £ A T implies that £\ A TC, £2 A TL 
or £i t 2 A H hold. Unfortunately, from (i29j) . (f3~Tj) and (|32l) we cannot guarantee that 
P(£ A J 7 ) is smaller than P(£ AT), but in any case, by multiplying these probabilities 
by [n]i in view of (f27l) . we complete the proof of (f24l) . 

Case 3 (s = w(l/(rn)) but also s = 0(r)). Following the same notation as in the 
case s = 0(l/(rn)) and by an analogous argument we obtain 

P(£A.F)~ (J±y (l- e -9«)4+4 e -^n_ ^y g -^n (33) 

If £3 < 1, we claim that this is the main contribution to P(£). In fact, suppose that TC 
holds and also that p > of the vertices in J are restricted (i.e. T does not hold). This 
happens with probability 0(r 2p ). Since £3 < 1, then the only possible event which 
contributes to S required in the definition of £ is (Si), involving vertex £ which cannot 
be restricted by definition. Then we deduce that Vol iJZ) > (I — p)irr 2 + £%q + eirr 2 , 
since the unrestricted vertices in J contribute {£ — p)irr 2 + £%q to Vo\(JZ) and the first 
restricted one gives the term enr 2 , by Lemma 13.21 (2—3) . Therefore, the probability of 
£ in this situation is 0(e~^ 3qn /n e -~ p+e ), which combined with the probability 0(r 2p ) 
that p vertices are restricted has negligible weight compared to fj33[) . Hence, P(£) ~ 
P(£ A J 7 ), and the first line of §2B follows from (l27j) and (|33l) . 

Unfortunately, if £3 = 2 and we have p restricted vertices in J, we can only assure 
that Vol (1Z) > (£ — p)irr 2 + q + eirr 2 , and then for some < e' < e 

P (£ A T) = O e~*> = O (--L) e-*\ (34) 

which may have significant contribution to P(£) if s is large enough. But in any case, 
in view of (|27|) . ([33]) and ([3"3]) . we verify that the second line of (f25j) is satisfied. 
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Case 4 (s = u{r)). Let T 1 be the event that for any i,j £ J (i ^ j) we have 
> 2r and also d(J^,X<) > 2r. This event has probability 1 - 0(r 2 ). We 
observe that if T' holds, then for any i,j £ J (i / j) we must have TZi fl 7£j = 0, 
ft', n 7^- = and Hi n 7^ = 0. Therefore, Vol (ft) = £vrr 2 + £ 3 q and the sets in Q are 

pairwise disjoint and also disjoint from ft. Then in view of Lemmata 13,31 and 13.41 and 
by the same argument that leads to (|28p 

P(S A T') ~ (1 - e -9")^l+^ e -4?n _ e -£ 3 <?n (35) 

The remaining of the argument is analogous to the previous case but replacing T with 
J- 1 and using Lemma 13. 2( 4). □ 

Taking into account that Kit = -D* + St and Kit+i = St + Bt, the number of 
isolated vertices at two consecutive steps can in the case s = 0(l/(rn)) be completely 
characterized by Proposition 13.11 For the other ranges of s, the result is weaker but 
still sufficient for our further purposes. We remark that if s = o(l/(rn)) then creations 
and destructions of isolated vertices are rare, but a Poisson number of isolated vertices 
is present at both consecutive steps. Otherwise if s = w(l/(rn)) then the isolated 
vertices which are present at both consecutive steps are rare since, but a Poisson 
number of them is created and also a Poisson number destroyed. 

Now in order to characterize the connectivity of (G(Xt\ r )) tg ^> we nee d to bound 
the probability that components other than isolated vertices and the giant one appear 
at some step. We know by Theorem 12.11 that a.a.s. this does not occur at one single 
step t. However during long periods of time this event could affect the connectivity 
and must be considered. 

Extending the notation in Section [21 given a step t let K%,t be the number of non- 
solitary components other than isolated vertices occurring at step t. We show that 
they have a negligible effect compared to isolated vertices in the dynamic evolution 
of connectivity. 

Lemma 3.5. Assume that jjL = 0(1) and s = o(l/(rn)) . Then, 

• P(K 2>t > A K 2tt+ i = 0) = P(K 2 , t = A K 2tt +i > 0) = o{srn), 

• P(K 2 , t > A B t > 0) = o(srn). 

Proof. Recall from Lemma 13.31 that if s = o(l/(rn)) then q = 0(rs). Then it is 
enough to prove that P(K 2 j > A K 2>t+ x = 0) = o(qn) and P(K 2>t > A B t > 
0) = o(qn), since (K 2> t = A K 2) t+i > 0) corresponds in the time-reversed process to 
(K 2> t > A K 2j t+i = 0) and thus they have the same probability. 

Consider all the possible components in G{X; r) which are not solitary and have 
size at least 2. They are classified into several types according to their size and 
diameter, and we deal with each type separately. Then if we denote by Mj the 
number of components of type i in G(Xt;r), we must show for each i that 

P(Mj > A k %t+ x = 0) = o{qn) and P(M t > A B t > 0) = o{qn). (36) 
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Also we need one definition which helps to describe the changes of edges between 
G(X t ; r) and G(X t+1 ; r). For each % G {1, . . . , n} we define Vi = Qi U Q- = 7^A7^ 
(where A denotes the symmetric difference of sets). Given also j G {l,...,n}, see 
that X,- € Vi iff Xj G "Pj iff vertices i and j share an edge either at time t or at time 
t + 1 but not at both times (which happens with probability Vol ("Pi) = 2q). 

Each part in this proof is labelled by a number followed by a prime (') in order 
to avoid confusion with the parts in the proof of Lemma 12.21 which are often referred 
to. Moreover, we write for simplicity Part i (p.L. I2.2j) to denote Part i in the proof 
of Lemma 12.21 

We set throughout this proof e = 10~ 18 . 

Part 1' . Consider all the possible components in G(X;r) which have diameter at 
most er and size between 2 and log re/37. Call them components of type 1, and let 
Mi denote their number at time t. This definition is similar to the one in Part 1 
(p.L. 12.2]) . but also includes components of size 2, covered by Lemma 12. 11 

Given any % G {1, . . . , re}, let £j be the following event: There exists a component 
r of type 1 in G(X \ {Xj}; r) and moreover for some j G {1, . . . , re} such that Xj is a 
vertex of T we have that Xj G Vj . In order to compute the probability of £ j , we note 
that the arguments in the proofs of Lemmata 12.11 and 12.21 are still valid if we replace 
X by X \ {Xj} (i.e. we ignore vertex i in the model). Hence, the probability of having 
some component in G(X \ {Xj}; r) of type 1 and size at least £ > 2 is 0(1/ log^ -1 re). 
Suppose first that G(X \ {Xj};r) has some component T of type 1 and size between 
3 and logre/37. This happens with probability 0(l/log 2 re). Conditional upon this, 
the probability that Xj G Vj for some j G {1, . . . , re} with Xj being a vertex of T is at 
most logre/37 times 2q. This contributes O (1/ log 2 re) (log re/37) (2q) = 0(q/ log n) to 
the probability of £j. Otherwise suppose that G(X \ {Xj}; r) has some component T 
of type 1 and size exactly 2. This happens with probability 0(1/ log re). Conditional 
upon this, the probability that Xj G Vj for some j G {1, . . . , re} with Xj being a vertex 
of T is at most two times 2q. This also contributes 0(1/ log re) (4q) = 0(q/ logn) to 
the probability of £j, and therefore P(£j) = 0(g/logre). 

Given any ii,i2 G {l,...,re} (i\ / 12), let be the following event: There 

exists a component T of type 1 in G(X\{Xi 2 }; r) and moreover TZ' ii n(X\{Xi 1 , Xj 2 }) = 
0. To derive the probability of J-i ± .j 2 , we distinguish two cases according to the distance 
between Xj 1 and V. Suppose first that for some h G {1, . . . , re} \ {ii, 12} we have that 
r < c^Xj^X/J < 3r (which happens with probability 0(r 2 ) = 0(log re/re)). Let Sh be 
the set of points in the torus [0, l) 2 at distance greater than er but at most r from X^, 
and let 5j x be the circle with center Xj x and radius r — 2s. At least one halfcircle of 
5j 1 has all points at distance greater than r from X^, so Area(5/ l U5j 1 ) > (1 — e 2 )irr' 2 + 
vr(r - 2s) 2 /2 > (5/4)7rr 2 . Notice that, if J zr j li j 2 holds for some component T which 
contains a vertex X^ such that d(Xj 1 ,X/ l ) < 3r, then we must have d(Xj 1; X/i) > r 
and moreover U must contain no point in X \ {Xj 1 ,Xj 2 }, which occurs with 
probability (1 — Area^ U cSjJ) 71-2 = 0(l/n 5 / 4 ). Therefore, multiplying this by the 
probability that (i(Xj 1 ,X/ t ) < 3r and also taking the union bound over the re — 2 
possible choices of h, the contribution to P(J 7 j li j 2 ) due to situations of this type 
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is O (n (log n/n)/n 5 / 4 ) = 0(log n/n 5 / 4 ). However, we claim that this has negligible 
asymptotic weight in P(J r j lj j 2 ). In fact, the probability that J-"i lt i 2 holds for some 
component T with all vertices at distance greater than 3r from is 0(l/(nlogn)) 
(in fact we only need an upper bound and therefore we just show that this probability 
is 0(l/(nlogn)). In order to prove this last claim, we consider all the notation in the 
proof of Lemma |2. II for the remaining of the paragraph, and also define S = tt^ 1 (S) 
and y = 7r- 1 (y). Then we can repeat the same computations there but, instead of 
asking that all the n — i points in X \ y lie outside of S, we require that all the 
n — i — 2 points in X \ (y U {X^ , Xi 2 }) lie outside of S U 1Z\ . This last fact occurs 
with probability P = (1 — Vol(5 U T*^))™ - ^ -2 , which plays a role analogous to that 
of P. If X{ x is at distance greater than 3r from any point in y, then S and TZ^ are 
disjoint. Therefore from (j3]) we get 



irr 



^(2 + ^)<Vol(5uK)<^ (37) 
and an argument analogous to that leading to © shows that 

P < (n) (l-137rr 2 /4)* +1 " (38) 

Then, repeating the same computations in the proof of Lemma 12.14 but replacing P 
with P, proves the claim for components of type 1 of fixed size I > 2, and this is 
extended to all components of type 1 by arguing as in Part 1 (p.L. 12. 2p . As a result, 
we conclude that P(J r j 1) j 2 ) = 0(1/ (n log n)). 

Now we proceed to prove (j36|) for components of type 1. First observe that the 
event (M\ > A i^2t+i = 0) implies that £,{ holds for some i £ {1, . . . ,n}, since the 
only way for a component of type 1 to disappear within one time step is getting joined 
to something else. Therefore, 

P(Mi > A K w = 0) < V P(Si) = O (t^~) ■ 

Also notice that (M\ > OABt > 0) implies that J~i lt i 2 holds and moreover X{ 2 G , 
for some i%, ii G {1, . . . ,n} (i\ ^ 12)- Then, 

P(M t > OA B, > 0) < EPft. A (X„ e Si,)) = O (^) =<>(£). 
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Part 2'. Consider all the possible components in G(X; r) which have diameter at most 
er and size greater than logn/37. Call them components of type 2, and let M2 denote 
their number at time t. 

Repeat the same tessellation of [0, l) 2 into cells as in Part 2 (p.L. 12. 2p . and also 
consider the set of square boxes defined there. Given any box b and i,j G {1, . . . , n} 
(i 7^ j), we define £b,i,j to be the event that box b contains more than logn/37 — 1 
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points of X \ {Xi} and moreover Xi 6 Pj. Observe that each of the events {M% > 
A A^2,t+i = 0) and (M% > A B t > 0) implies that £b,i,j holds for some box b and 
i,j £ {1, . . . , n}. Then, by repeating the argument in Part 2 (p.L. 12.2]) . but ignoring 
Xi and also replacing logn/37 with logn/37 — 1, we deduce that 

P(M 2 > A K w = 0) < O (-^- ) V P(Xj 6P») = f off ) , 

\n LI logn/ ^— ' \n Ui logn/ 

and the same bound applies to P(M2 > A fij > 0). 

Pari 3'. Consider all the possible components in G(X;r) which are not embeddable 
but not solitary either. Call them components of type 4, and let M4 denote their 
number at time t. 

Repeat the same tessellation of JO, l) 2 into cells as in Part 4 (p.L. I2.2h . and observe 
that each of the events (M 4 > A K 2) t+i = 0) and (M 4 > A B t > 0) implies that for 
some i,j £ {1, . . . , n} there exists some connected union S* of cells in the tessellation 
with Area(S*) > (ll/5)vrr 2 such that S* n (X \ {Xi}) = and moreover Xi G Vj. 
Hence, from Part 4 (p.L. 12.2]) but replacing X with X \ {Xi}, we obtain 

P(M 4 > A K W = 0) < O (-^ ) p (^' = ( J™ ) , 

\n°'° log nj ~ \n L ' logn/ 

and the same bound applies to P(M 4 > A B t > 0). 

Part 4'. The embeddable components with diameter at least er treated in Part 3 
(p.L. 12. 2p are here divided into two types. First consider all the possible components 
in G(X;r) of diameter between er and 6\A2V- Call them components of type 3a, and 
let M^a denote their number at time t. 

We tessellate the torus [0, l) 2 into square cells of side ar, for some fixed but 
small enough a > 0. From Part 3 (p.L. [272]) . if G{Xt\r) has some component of this 
type, then there exists a topologically connected union S* of cells with Area (S*) > 
(1 + e/6)7rr 2 which contains no point in X. By removing some extra cells from 5*, 
we can assume that the number of cells in S* is exactly |" ^ 1+ ^ 6 ^ 7r ] . Now for each 
i, j £ {1, . . . , n} and each union S* of [" ( 1 + e / 6 ) 7r "j ce \\ s which is topologically connected, 
let £i,j,s* be the following event: S* contains no points in X \ {Xi,Xj}, Xj is at 
distance at least 2r from all the points in S*; TZj contains no points in X \ {Xi,Xj}; 
and moreover Xi £ Vj. Notice that if Xj is at distance at least 2r from all the points 
in S*, then ttT" 1 ^*) and VJ j are disjoint. Hence, Vol^T" 1 (S*) U H^) > (2 + e/6)vrr 2 
and 

P(^-,5*) < (l - Vol^OS*) U TZ^y 2 (2a) = O (-^) 

Similarly, let J~i,j,s* be the following event: S* contains no points in X \ {Xi,Xj}; 
Xj is at distance at most 2r from some point in S*; and moreover Xi £ Vj. Notice 
that the probability that Xj is at distance at most 2r from some point in S* is 
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0(r 2 ) = 0(logn/n). Hence, 



P<^.) < (1 - mr ->0 (^) ,2„ = o (^J) 

Finally, observe that each of the events (M 3a > A K 2 ,t+\ = 0) and (M 3a > A B t > 
0) implies that either £i.j,s* or J-"i,j,s* hold, for some i,j G {l,...,n} and some 
topologically connected union S* of cells. Therefore, the probabilities of (M 3(l > 
A K 2 , t +i = 0) and (M 3a >0A5 ( >0) are at most 

Pari 5'. Finally consider all the possible components in G(X; r) which are embeddable 
and have diameter at least 6\/2r. Call them components of type 3b, and let M 3 ?, 
denote their number at time t. 

We tessellate the torus into square cells of side ar, for some fixed but small enough 
a > 0. Our goal is to show that if G(Xt; r) has some component of type 3b, then there 
exists some topologically connected union S* of cells with Area(5*) > (ll/5)7rr 2 and 
which does not contain any vertex in X . Then, arguing as in Part 3', we conclude 
that both P(M 36 > A K 2 ,t+i = 0) and P(M 36 > A B t > 0) are O (qn/in 1 / 5 log n)) . 
We now proceed to prove the claim on the union of cells S* . Given a component V 
of type 3b in G(X t ;r), let 5', ij and ie be defined as in Part 3 (p.L. I2.2[) . Then, by 
repeating the same argument in there (but replacing er with 6s/2r), we can assume 
wd.o.g. that the vertical distance between Xi T and Xi B is at least 6r, and claim 
that the upper halfcircle with center Xi T and the lower halfcircle with center Xi B 
must be disjoint and contained in S' . Now, consider the region of points in the torus 
[0, l) 2 with the y-coordinate between that of Xi T and Xi B , and split this region into 
three horizontal bands of the same width. Observe that each band has width at 
least 2r and hence must contain some vertex of T. For each of these bands, pick the 
rightmost vertex of T in the band. We select the right lower quartercircle of radius 
r centered at the vertex if the vertex is closer to the top of the band, or the right 
upper quartercircle otherwise. We also perform the symmetric operation and choose 
three more quartercircles to the left of the leftmost vertices in the three bands. All 
this six quartercircles together with the two halfcircles previously described are by 
construction mutually disjoint and contained in S'. Therefore Area (5') > (5/2)-7rr 2 . 
Let S* be the union of all the cells in the tessellation which are fully contained in 
S' . We loose a bit of area compared to S' . However, if a was chosen small enough, 
we can guarantee that S* is topologically connected and also Area(5*) > (ll/5)-7rr 2 . 
This a can be chosen to be the same for all components of type 3b. □ 

Now we can characterize the connectivity of [G(Xt\ r)) t z at two consecutive steps. 
We denote by Ct the event that G{X t ; r) is connected, and by V t = Ct the event that 
G(X t ;r) is disconnected. 
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Corollary 3.1. Assume that fj, = 0(1). Then, 

P(C t A V t+1 ) ~ e -A«(l - e - EB ), P(p t A C m ) ~ e^l - e~ EB ) 
P(C t A Ct+i) ~ e-^- EB , P(P t A P m ) ~ 1 - 2e~^ + e~^ B 

Proof. First observe that K\ t t = St + D t and = St + Sj. Therefore we have 

P{K U = A > 0) = P(5 t = A A = A B t > 0), 

and by Proposition 13.11 we get 

P(K lft = A K w > 0) ~ e- E5 - ED (l - e- EB ) ~ e-"(l - e- EB ). (39) 

We want to connect this probability with P(C t A V t+ i). In fact, by partitioning 
[K\ t t = A Kij+i > 0) and (Ct A V t +i) into disjoint events, we obtain 

P(K ljt = A A\ m > 0) = P(C* A K w > 0) + P(P t A K u = A # 1)f+1 > 0), 
P(Ct A A+i) = P(Ct A ifi )t+ i > 0) + P(Cj A 2? t+ i A JRTi,*+i = 0), 

and thus we can write 

P(Ct A A+i) = P(#i,t = A A w > 0) + Pi - P 2 , (40) 

where P x = P{C t A V t+1 A A\ m = 0) and P 2 = P(D t A K u = A A\ t+1 > 0). 

Now suppose that s = o(l/(rn)). In that case, P(Ki jt = OAi^i^+i > 0) = Q(srn) 
(see ([39]) and Proposition 13 . X [> . Also observe that V A (X = 0) implies that X > 0. 
In fact, we must have at least two components of size greater than 1, so at least one 
of these must be non-solitary. Then, we have that Pi < P(K 2 t = OA^^i > 0) and 
P 2 < P(A\t >0A5(>0), and from Lemma [331 we get 

Pi,P 2 = o(P(Ai ii = 0AAi, m >0)). (41) 

Otherwise if s = fi(l/(m)), then P(Ai,t = A Ai )t+ i > 0) = 0(1). In this case, we 
simply use the fact that Pi < P(K 2>t+ i > 0) = o(l) and P 2 < P(K 2jt > 0) = o(l) 
(see Theorem 12.21 and Lemma [3 . 1 f) . and deduce that (141D also holds. 

Finally, the asymptotic expression of P(Ct A 2?t+i) is obtained from ([39]) . (|4"0"]) 
and (f4T) . Moreover, by considering the time-reversed process, we deduce that P(V t A 
Ct+i) = P(Cf APt+i). The remaing probabilities in the statement are computed from 
Corollary 12.11 and Lemma 13.11 and using the fact that 

P{C t A C m ) = P(C t ) - P{C t A V t+1 ), 

P(V t AT> t+1 ) = P(p t )-P(V t AC t+1 ). □ 

Let A be an event in the static model G(n,r). We denote by At the event that 
A holds at time t. In the [G(Xt]r)) t z model, we define Lt(«4) to be the number 
of consecutive steps that A holds starting at step t (possibly if A t does not hold). 
Note that the distribution of L t (A) does not depend on t, and we will often omit the 
t when it is understood or not relevant. 
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Lemma 3.6. Consider any event A in the static model. If we have that E(L(A)) < 
+00 (but possibly E(L(A)) — ► +00 as n — ► i/ien conditional upon At but not 

A t -\ we have 

E(L t (A) I A^i A At) = vr ^l tt , 

P(A-i a A) 

which does not depend on t. 
Proof. We have that 

Lt_i + l|AIT]L t = l[A-i] + U 
and taking expectations and using the hypothesis that E(L(«4)) < +00 we get 



B(l[A t -i]L t (A)) = P(A), Vt. 

Using the fact that 



E(L t (A I A-l A A t 



E(l[A-i A A]MA) E(l[A t -i]L t (A)) 



P(A-i a A) P(A-i a A) 

the result follows. □ 

To prove that E(L(C)) < +00 and E(L(£>)) < +00 we need the following technical 
lemma. 

Lemma 3.7. Let b = b(n) be the smallest natural number such that (b — 3)ms > 
3v2/2. Then, there exists p = p(n) > sitc/i t/iai: for any fixed circle 1Z C [0, l) 2 
of radius r/2, any i € {1, . . . ,n}, any t G Z, and conditional upon any particular 
position of X, Lt t in the torus, the probability that X^t+bm is at least p. 

Proof. Fix an arbitrary position for Xij and also for circle 1Z, and call X to its center. 
Let t' be the smallest integer such that t' \ m and t' > t (i.e. t' is the first time after 
t when agent i selects a new angle), and call h = t' — t, which naturally satisfies 
< h < m. The particular point X^t is irrelevant in our argument, and we restrict 
our attention to the position of agent i at the times when it chooses a new angle 
(plus the final position), and call for simplicity = Xi^+km (0 < k < b — 1) and 
Y;, = Xij+bm. Observe that 

d(Y k+1 , Y k ) = ms, Vk : < k < b - 2, and d(Y b , Y^) = (m - h)s. (42) 

Moreover recall that, if a& denotes the angle in which agent i moves between Y k 
and then each a k is selected uniformly and independently at random from the 

interval [0, 2ir). 

In order to prove the statement, we compute a lower bound on the probability of 
a strategy that is sufficient for agent i to reach 1Z at time t + bm. We start from an 
arbitrary point Yq G [0, l) 2 and build a sequence of points Yq, satisfying (|42|) 

such that d(Yb,X) < r/2, by imposing some restrictions on the angles cto, ...,a%. 
For the sake of simplicity in the geometrical descriptions, it is convenient to allow 
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Yq, . . . ,Yb and X to lie in M? rather than into the torus [0, l) 2 . Once the construction 

of the sequence of points is completed, we map them back to the torus by the usual 

projection. Hence, we assume hereinafter that Yq and X are two arbitrary points in 

M. 2 such that d(Yo,X) < v2/2 (which is the maximal distance in the torus [0, l) 2 ). 

For each k, < k < b — 4, we restrict a k to be in [Ok — 7r/6, Ok + vr/6] (mod 2ir), where 
> 

Ok is the angle of Y k X with respect to the horizontal axis. We claim that, with this 
choice of angle, the distance between Yk and X is decreased at each step by at least 
ms/3 until it is at most ms. In fact by the law of cosines, 



d(Y k+1 ,X) < y {d(Y k ,X)) 2 + (ms) 2 - Vsd(Y k ,X)ms, (43) 
and therefore, if d(Y k , X) > ms, we can write 



d(Y k+1 ,X) < x j (d{Yk,X)f + (l + ~ Vt)(ms) 2 - ^d{Y k ,X)ms 



< \j (d(Y k ,X)) 2 + ^(ms) 2 -^d(Y k ,X)ms 



= d(Y k ,X)-~ms. (44) 
Otherwise, if d(Y k ,X) < ms, then from (|4"3|) we deduce that also 



d(Y k+1 , X) < yJ(l-V3)(d(Y k ,X)) 2 + (ms) 2 < ms. (45) 

Hence, we can guarantee that d(Yb_3, X) < ms: Suppose otherwise that d(Y&_3,X) > 
ms. Then in view of (|43[) . (|44|) and (|45p . for all k such that 0<fc<6-4we also 
have d(Yk,X) > ms, and moreover 

d{Y b ^,X) < d(Y ,X) - (6 - 3)^ < ^ - (6 - 3)^ < 0, 

which contradicts the assumption. 

Let Z S M. 2 be the only point on the line containing Y5-3 and X satisfying 
d(Z, X) = (m — h)s and such that X lies on Y^-^Z. Denote by W one of the two 
points on the perpendicular bisector of Yt,_ 3 Z which satisfy d(W, Y&_3) = ms. We 
want to set the angles a&_3, aj,-2 and so that 1&_2) ^fe-i and ^6 are close to 
W, Z and X, respectively. Indeed, if ^-3, 4>b-2 and 4>b-i, respectively, are the an- 
gles between the horizontal axis and Y^-^W , WZ and ZX, then by imposing that 
o k £ [4>k — er/ (ms), 4>k + er/(ms)] (mod 2ir) for some small enough e > 0, we achieve 
that d(Y h ,X) < r/2 and thus Y b G TZ. 

Therefore, the probability of choosing all the angles according to the strategy 
described is p = (l/6) 6 ~ 3 © ((r/(ms)) 3 ). □ 

The next lemma allows us to apply Lemma 13.61 

Lemma 3.8. E(L(C)) < +00 and ~E(L(T>)) < +00. 
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Proof. Fix one circle TZ C [0, l) 2 of radius r/2, and take b as in the statement of 
Lemma 13.71 Since the agents choose their angles independently from each other and 
in view of Lemma f3.71 we have that, conditional upon any arbitrary Xt, the probability 
that all agents end up inside TZ after bm steps is 

P(X t+bm <zK\Xt)> p n , (46) 

for some p = p(n) > 0. Observe that for any f £ Z the event (Xt C TZ) implies that 
G(X t ; r) is a clique, since all pairs of vertices in X t are at distance at most r, and thus 
G(X t ; r) is connected. Consequently, for any d E N, we can write 

d-l 

n\d 



( A Vt+kbm) < (1 "P n )p( A Vt+kbm) < P(A)(1 
^fc=0 ' ^ k=0 '' 



p n f. (47) 



Now observe that the equation L t (V) = X^fcLo l[£>t] • • • lp^t+fc], is satisfied point- 
wise, for every element in the probability space (Xt)t&- Therefore, by the Monotone 
Convergence Theorem, (|47j) and the fact that p > 0, we conclude 

oo 

E(L t (2?)) = ^P(P t A---A2? m ) 

fc=0 

oo / d \ 

+kbm I 

d=0 ^ A:=0 ' 

oo 

< bmV(V t )^(\ - p n ) d < +oo. 



d=0 



The same kind of argument applies to show that E(L(C)) < +oo. In this case 
we fix two circles TZ and 1Z' in [0, l) 2 of radius r/2 and far apart from each other 
(say with centers at distance greater than 2r). Observe that for any t G Z the event 
((X t \ {Xij}) C TZ) A (Xi tt G 7£') implies that G(^;r) is disconnected. Moreover, 
from Lemma 13.7} we obtain an analogue to (1461) 



\ {X u+hm }) C 72.) A (X 1)t+6m G 72') | > p», (48) 

and the argument follows as in the previous case but replacing T> with C. □ 

We are now ready to prove our main theorem which characterizes the expected 
number of steps the graph remains (dis) connected once it becomes (dis)connected. 

Theorem 3.1. 



B(L k (C) | V k ^AC k 



(1-e 



ifsrn = o(l), 
(i- e -LnA) ifsrn = e(l), 
1 if srn = u(l), 
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( Tr(e^-l) 



E(L k (V) | C fc _i AP fc ) 



1 



(1 



Asm 



(1-e- 



i/ srn = o(l), 
-w^) ifsrn = Q(l), 
1 i/ srn = w(l). 



Proof. Since by Lemma [3,81 E(L&(C)) < +00, E(Lfc(P)) < +00, we can apply the 
formula of Lemma 13.61 and the results follow by Corollary 13.11 □ 



4 Conclusion. 

In this paper we have formally introduced the dynamic random geometric graph in 
order to study analytically dynamic MANETs. We studied the expected length of 
the connectivity and disconnectivity periods, taking into account different step sizes s 
and different lengths m during which the angle remains invariant, always considering 
the static connectivity threshold r = r c . We believe that a similar analysis can be 
performed for other values of r as well. 

The Random Walk model simulates the behavior of a swarm of mobile vertices as 
sensors or robots, which move randomly to monitor an unknown territory or to search 
in it. There exist other models such as the Random Way-point model, where each 
vertex chooses randomly a fixed way-point (from a set of pre-determined way-points) 
and moves there, and when it arrives it chooses another and moves there, and so on 
[3] . A possible line of future research is to do a study similar to the one developed in 
this paper for this way-point model. We believe that the techniques developed in this 
paper will prove very useful to carry out that study. 
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